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We present a general theory to describe two-photon interference, including a formal description of 
few photon intereference in terms of single-photon amplitudes. With this formalism, it is possible to 
describe both frequency entangled and separable two-photon interference in terms of single-photon 
wave functions. Using this description, we address issues related to the physical interpretation of 
two-photon interference experiments. We include a discussion on how few-photon interference can 
be interpreted as a bosonic exchange effect, and how this relates to traditional exchange effects with 
fermions. 

PACS numbers: 42.50.Dv, 03.65.-w 

Newton and WignerQ first discussed the non-localizability of photons, which prevents the introduction of a position- 
representation wave function in the usual sense of a wave function for a massive particle. Nonetheless, it has been 
known for years that it is possible to present physically meaningful descriptions of photon detection in finite regions 
of space 2] . This has led to a host of approaches to introduce a "photon wave function" Q , each of which establishes 
at least a limited analogy to massive particle wave functions. Today there is no single accepted definition of such 
a photon wave function, as there are a variety of possible analogies, and the most convenient often depends on the 
system one would like to describe. For example, Chan et al. discuss a wave function associated with a photon 
spontaneously emitted from an atom in terms Schmidt pairs of atomic and photonic cigcnfunctions, whereas Resch et 
al. p| find it useful to consider a photon wave function based on the Glauber detection probability in understanding 
an absorptive exchange effect. Given the variety of convenient types of "photon wave functions" introduced in the 
literature, it is perhaps best not to insist on a particular definition, but rather to understand the term to refer broadly 
to any approach for describing a photon in a manner analogous to the usual massive particle wavefunctions introduced 
in nonrelativistic physics. That is the point of view we take here, where we use photon wave functions to describe 
few-photon interference experiments. Although in the particulars of our discussion we use a photon wave function 
definition based on the Glauber detection formula, within the usually relevant approximations an easy translation 
into other photon wave function definitions could be made. 

Yet one might ask, "why bother?" After all, the measurement results of two-photon interference experiments can 
be predicted with relatively simple, well-known calculations. Still, what seems to be a less straightforward task 
is the physical interpretation of the experiments. Few-photon interference is often discussed loosely in terms of 
interfering Feynman paths, or the overlap of wave packets, or the distinguishability of a particular set of outcomes. 
In this paper we formally address what one must mean if one wishes to discuss the interference of individual photon 
amplitudes leading to a particular detection event. Using photon wave functions, we explicitly show that second-order 
interference experiments can be understood in terms of single-photon amplitudes. Consequently, this approach yields 
considerable insight into few-photon interference experiments, and illustrates the relevance of the bosonic nature of 
the photon[3| by contrasting the corresponding description that would result for fermions. It is with a photon wave 
function description that one can perhaps best isolate non-classical interference terms, and discuss the interference of 
single-photon amplitudes corresponding to a particular detection. 

Perhaps the most familiar few-photon interference experiment involves the Hong-Ou-Mandel interferometer We 
show how second-order interference of this type manifests itself as an exchange effect in the photon wave function 
picture. This seems to be appreciated by many workers in the field, although we have not been able to find an explicit 
discussion of this point in the literature. We then move on to other experiments. We show how the elimination of 
which-path distinguishing information restores the exchange effect in quantum eraser experiments |l l| . And while in 
their discussion of a postponed compensation experiment |12| Pittman et al. emphasize the limitations of using a pair 
of photon wave packets to describe frequency entangled two-photon interference, we show that it is indeed possible to 
understand this experiment in terms of the interference of single photon amplitudes. Before proceeding to these issues 
in section II below, we identify our definition and notation for few-photon wave functions in section I; our conclusions 
are presented in section III. 
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I. PHOTON WAVE FUNCTIONS 



Photon wave functions can be extracted from the usual field theory description of photodetection, and used in a 
way similar to the use of massive particle wave functions in describing the detection of those particles. To establish 
this analogy we recall some well-known results from the field theory for nonrelativistic massive particles, in which an 
arbitrary single-particle state \S) can be written as 

\S)=J29 S W^W\vac), (1) 

A 

where each A identifies one of a set of normalized, orthogonalizcd single-particle modes, \vac) is the vacuum state, 
and at (A) is the associated particle creation operator. For particles in free space we can take A = (//, k) to label both 
the spin state \i and the wave vector k of a plane wave. Here g (A) is a normalized amplitude, satisfying 

$> S (A)f=l. (2) 

A 

We can write this state JJJ in terms of the creation operator for a particle with spin label /i at r, 



i e lk r 



where V is a normalization volume, as 



where 



|5)=E/ drGS(r)*l(r)\vac), (3) 



G^(r) = JL]r/(AK k - r , 

and we anticipate the passage to the limit of infinite V and the resulting continuous range of k over plane wave 
modes. We can now identify the single-particle wave function <5?(r, t) associated with the single particle state 15*). 

The state $ M (r,i) IS"), where $ M (r,t) = U (i, 0) $ M (r) £/t (t, 0) and U (t, 0) is the time evolution operator from to 
t determined by the Schrodinger equation, is proportional to the vacuum state; the wave usual function $^(r,t) of 
elementary quantum mechanics provides just that proportionality, 

$ p (r,t) \S) = (r,t) \vac) , (4) 

as can be easily confirmed. The probability of detecting the particle within dr of r at time t can now be written in 
either the particle or field theory notation, 

]T (r,t)f dr = (S\ *], (r,t) $ M (r,t) |5) dr. (5) 
Moving to two-particle states, we can construct the most general such state \T) according to 

\T) =N T J2 f A') at (A) at (A') \vac) = E / / d ^' F l' ( r - r ') H ( r ) % ( r ') l« ac ) » ( 6 ) 

A, A' MiM 7 

where for later convenience we have introduced an explicit normalization constant iV T such that (T\T) — 1, and 

<^(r,r') = ^E/ T ( A ' A ')e ik V k ' r '. 

k.k' 
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Without loss of generality the amplitude f T (A, A') can be taken to be symmetric with respect to the interchange of 
A and A' for bosons, and antisymmetric with respect to that interchange for fermions. A two-time, two-particle wave 
function (or bi-particle wave function), <J>^ , (r, r',t, t'), can be introduced according to 



-j=% (r,t) V (r',0 \T) = ^ (r, r', t, t') \vac) . (7) 

It is easy to show that ^„„/ (r, r', t, t) is the usual two-particle wave function at time t, symmetric (or antisymmetric) 
with the interchange of (r,ju) and (r', /i'), if we deal with bosons (or fermions). And considering detection processess 
activated at time t, a standard calculation easily done at the wave function level shows that the probability of detecting 
one particle within dr of r and a second particle within dr' of r' is given by 

(r,r', M)| 2 drdr' = ^ £ ( T % (r',i) % (r,t) $ M (r,t) $ M , (r',t) r) drdr'. (8) 
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We now introduce photon wave functions in such a way that the equivalences (|5I8|) between wave functions and field 
theory descriptions hold when the standard Glauber detection formulas are used to model photodetection probabilities. 
Between detector activations we assume that the electromagnetic field evolves as a free radiation field. A modified 
version of this approach can be written down if this is not the case, but we will not do so here. We begin by introducing 
a single-photon state by analogy with (JJJ , 

|S> = 5> s (A)at(A)|mc) (9) 

A 

where (A) is the photon creation operator for a mode A, where the index A labels both the polarization (or helicity) 
and the wave vector k, and the amplitudes {g s (X)} again satisfy (J2J. The single-photon wave function at position r 
and time t, <& s (r,t) is defined by 

E+(r,t) \S) = * s (r,i) | vac) (10) 

where 

E+(r,t) = z£x fc e A a(A)e l ( k ™ t ), (11) 

A 

E"(r,t) = (E+frt)) 1 

are respectively the positive and negative frequency components of the electric field operator E(r, t) = E + (r,i) + 
E~(r, t), Kk = \j2i:hujk/V with V as the normalization volume, e\ is the unit vector indicating the polarization, and 
Uk = c |k| is the angular frequency of wave vector k. The analogy of this photon wave function with the single particle 
wave function of a nonrelativistic particle lies in the fact that, for the single photon state, the Glauber detection rate 
of a single ideal detector at position r and at time t is proportional to the first order correlation function, given byQ 

w^(r,t) = (5|E-(r,t)-E+(r,t) \S) = |* S (r,i)f, (12) 

cf. (Sj. For photon wave functions this is, of course, not a fundamental postulate of the theory, but is derived from 
the Glauber detection model; in particular, we use a form of the model where the detector is assumed isotropic. The 
dynamics of $ s '(r, t) are not governed by the Schrodinger equation, as in the massive particle case, but rather by 
Maxwell's equations. 

The analogy with massive particles carries on to states with more than one excitation. A general two-photon state 
of the radiation field is described by 



\T) =N T Y,f (A, A') at (A) a t (A') \vac) (13) 

A, A' 

where N T is the normalization factor, which for photons (or massive bosons) can be taken to be 

^ T = (El/ T ( A ' A ')r+E[/ T ( A '' A )]*/ T (A,A')] , (14) 



,A,A' A, A' 
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to guarantee (T|T) = 1. We do not choose the function f T itself to be normalized, as it is convenient to write the 
normalization factor separately when comparing two-photon states to corresponding single-photon states. A special 
case is that for which the amplitude f T (A, A') in some basis (taken here to be that of polarizations or helicities and 
wave vectors) can be written as the symmetric product of a function of A and a function of A'. Each of these functions 
can then be associated with a one-photon state. In such a case we write / T (A, A') as the function f p ® (A, A') = 
g P {A)gQ (A') -I- g p (A')g® (A), where {g p (A)} and {g^(A)} are normalized spectral amplitudes J5J associated with 
single photon states, and we call the two-photon state separable. As is usually done in the literature, one can describe 
the same state with the simpler non-symmetric function f p Q (A, A') = g p (A)g® (A') since the additional component 
does not change the state (11311 . With this choice of f p Q (A, A'), for such separable states, 



N PQ 



£[<7 P (A)]V(A) 



(15) 



and N p Q = 1 only if the two single photon amplitudes are orthogonal; if P = Q we have a state of two identical 
photons, and N PQ = 

We introduce a two-photon wave function (sometimes called a biphoton^l) that satisfies 
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E+ (r 2 , t 2 )E+ (n , h ) \T) = $1 (n , r 2 , t u t 2 ) \vac) 



(16) 



where Roman subscripts denote Cartesian components, and $ T (ri,r2,ii,i2) are the components of the wave function 
associated with the state \T). In a general two-photon state, the coincidence detection rate of two ideal detectors at 
positions n, r 2 and at times ti, t% is proportional to the second order correlation function^ 



wW(r 1 ,r 2 ,t 1 ,t 2 ) = ^<r|Sr(r 1 ,t 1 ) J Br(r 2 ,i 2 ) J B+(r2,t2) J B+(ri,t 1 )|r) 

ij 

= 2^|$ p (r 1 ,r 2 ,ii,i 2 ) 2 , 



(17) 



cf. ©, where $£(ri, r 2 , ii, t 2 ) is symmetric under exchange of (i,ri,ti) with (j,r 2 ,i 2 ). 



II. INTERFERENCE EXPERIMENTS 



We now examine a type of quantum interference, observed in a number of few-photon experiments, that is associated 
with the measurement of particular coincidence detection rates. We discuss these rates in terms of interfering single- 
photon amplitudes, and show that the interference can be associated with photon exchange effects. 

The simplest type of two-photon state is a separable one, in which we can write f p Q (A, A') = g p {A)g® (A'). Here 
we can say that one photon has the spectral properties of 'Photon P'' and one photon has the properties of 'Photon 
Q' referring to the single-photon states \P) and \Q). For massive fermions, the analogous kind of separable spectral 
amplitude corresponds to a standard Hartree-Fock, single-determinant wave function. In many-body physics one 
usually characterizes such a state as free of the "correlation effects" that arise due to electron-electron interactions in 
more sophisticated models of the full, many-electron wave function. Nonetheless, there are dynamical consequences 
due to exchange effects even in separable states, which arise for bosons as well, and in particular photons. After 
considering the separable case, we generalize the discussion to address frequency entangled few-photon interference 
experiments that explicitly involve correlation effects. 

The two-photon wave function for a separable two-photon state takes the form 

N PQ r i 
(n, r 2 , t u t 2 ) = -7=- |*f (ri, h)$f (r 2l t 2 ) + if (r x , ti)#f (r 2l f 2 )j , (18) 

with the + arising instead of a — because photons are bosons, and where $f (r, t) and 3>?(r, t) are the single-particle 
wave functions associated with photons P and Q respectively. It is in this form that one can discuss interference in 
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terms of individual photon amplitudes. The second-order correlation function is 



w 



(2) 



(ri,r a , t^ta) = £ |iV PQ 



$f(ri,ti)| $ Q (r 2 ,i 2 ) + ^(rx,^) |*f(r a ,t 2 ) 

+ [ffCn.tx)]**?^!,*!) [$?(r 2 ,i 2 )] V(r 2 ,< 2 ) 
+ [^(ri.tx)]"^^!,*!) [$«(r a ,t a )l**«(r a ,f 3 ) . 



(19) 



The first two terms on the right-hand side of (|19fl are the classical independent-particle terms. If only these terms 
were present, the detection coincidence would be simply identified with the alternatives "photon P at (ri,ti) and 
photon Q at (r 2 , t a )" or "photon Q at (ri, t-y) and photon P at (r a , f a )" characteristic of independent detection events. 
The last two terms on the right-hand side are the interference or exchange terms. One can see that the exchange 
terms are proportional to the indistinguishability, or overlap, of the P and Q single-photon amplitudes in the region 
of interest. Non-zero exchange terms are indicative of the "indistinguishable Feynman paths" or "overlapping single- 
photon amplitudes" sometimes mentioned in the few-photon literature. 

It is important to consider the finite detection window of a realistic detector. Though the exchange terms in l|19|) 
may be non-zero for some values of t\ and t 2 , the exchange terms may integrate to zero over a finite detection time. 
Particularity if one were to consider photons in different frequency ranges, one would find that spatially overlapping 
single-photon amplitudes would give rise to instantaneously non-zero exchange terms that would integrate to zero 
over a realistic detection time. It would be possible to observe interference between photons in different frequency 
ranges provided that the detection window were sufficiently short, but in practice no interference would be observed 
regardless of the spatial configuration of the photon wavefunctions in the region of detection. In contrast, if the 
photons have the same peak frequency, the spatial overlap of single-photon amplitudes in the detection region is 
sufficient to give rise to exchange terms that do not integrate to zero over a realistic detection time. In the subsequent 
discussion we assume that we are dealing with photons that have the same peak frequency, and thus where there are 
overlapping single-photon amplitudes in the detection region there are exchange effects in coincidence measurements. 

We now apply this formalism to several important two-photon interference experiments, and demonstrate that in 
each case the interference can be understood as single-photon wave function amplitudes giving rise to an exchange 
effect. 



A. The Hong-Ou-Mandel interferometer 



First we consider the simple case of the Hong-Ou-Mandel interferometer. A schematic diagram of this experiment 
is shown in Fig. 1, where we label the input ports a and b, and the output ports c and d. Two photons are incident on 
the 50/50 beam splitter, with one of the photons delayed in time by St. We use a simplified notation, taking only one 
linear polarization associated with each port, and considering small wave packets with a narrow spread of wave vector 
components perpendicular to the direction identified by the port. Hence a one-dimensional treatment is possible for 
each port. We use, for example, a£ CO to denote a creation operator associated with a wave vector of magnitude k 
and in the propagation direction relevant for port b. In this shorthand the initial separable two-photon state, at a 
time just before the photon wave functions impinge on the beam-splitter, can be written as 



\P a , Q b )=J j dkdk'g p (k)gQ(k')al (k) a\ (k') \vac) , (20) 

where we have gone to a continuous range of wave numbers. In the full notation used in the previous section the 
functions g p (k) and g® (k) would of course be orthogonal, since the single-photon wave functions do not overlap. But 
in our shorthand notation we take photons P a and Qb to be the same when referenced to their own ports, with only 
the second delayed by a time St from the first; thus g®(k) = g p (k)e lu)kSt . In order to reduce the calculation to simple 
integration over one-dimensional coordinates we use the scalars, x n , and mode indices, j n , where x n is the spatial 
coordinate for the input port j n , together to identify a position r of interest. With this notation and Eq. (|2U|I the 
two-photon wave function in the region just before the beam splitter can be written as 



<S>(ji,x 1 ;j 2l x 2 ;t 1 ,t 2 ) = -= [^ P {ji,x 1 - 1 t 1 )^ Q (j 2 ,x 2 - 1 t2) + ^ Q (ji,x 1 ;t 1 )<P p (j 2 ,x 2 ;t 2 )] , 



(21) 
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where the single-photon wavefunction Q p (j, x;t) is defined by the obvious simplification of (|10fl and 

$ p (a,x;t) = J dkg p (k)e l ( kx ~^ =V p {x,t), (22) 

$ Q (b,x;t) = J dkg Q {k)e^ kx ^ kt ^ = V Q {x,t), 
$ p (6,x;t) = <f> Q (a,x;t) = 0. 



The origins of ports a and b are equally distant from the beam splitter, as shown in Fig. 1. We have given the one- 
dimensional single-photon wave functions in this input region a special label, V p ^ (x,t), so that we can illustrate 
how these input wave functions interfere in the detection region. The wave function amplitude corresponding to a 
two-photon detection event in the input region before the beam splitter is 

$(a,x 1 ;b,x 2 ;t 1 ,t 2 ) = x 2 ; a, Xl ; t 2 , t t ) = [V p (x u h)V Q {x 2 , t 2 )\ (23) 

= -±= [V P (xi,h)V p (x 2 ,t 2 -8t)] 
V2 

and, as one would expect, the square modulus of this function does not give rise to exchange terms and there is no 
interference in this region. We now wish to evaluate the two-photon detection rate in the output ports c and d; the 
origin of each of these ports is taken at the beam splitter. As is well known, the effect of the beam splitter is to effect 
a canonical transformation on the port operators^, 



4 (k) 



4 (k) 



"7T 



al(k)-iat(k) 
a\ (k) - ia\ (k) 



(24) 



where A is the time required for light to travel from the origins associated with ports (a, 6) to the beam-splitter, and 
it will be convenient below to use delayed time coordinates tx(2) = ^i(2) — A. The two-photon wave function in the 
region of the detectors has the following components 

$(c, xi; c, x 2 \ ti, t 2 ) = — ^= [$ P (c, xi; ii)$ Q (c, x 2 ;t 2 ) + $ Q (c, x x ; h)<i> p {c, x 2 ;t 2 )] 
2v 2 

$(d,xi;d,x 2 ;t 1 ,t 2 ) = — y= [® p (d, Xl ;ii)<f> Q (d,xr,t 2 ) + $ Q (d, xi;tx)<Z> p (d, x 2 ;t 2 )] 

<$>{c,xi\d,x 2 ;t x ,t 2 ) = <${d,x 2 ;c,xx;t 2 ,ti) = [$ P (c, xf, ii)$ Q (d, x 2 ; t 2 ) + $ Q (c, Xx] Zi)<I> P (d, x 2 ; t 2 )] . 

2v 2 

We can write the wave function in terms of the single-photon input amplitude since 

<t> p {c,x;t) = -iV p ( Xl t) 



(25) 



$ p (d,x;t) = V p {x,t) 

$ Q (c,x;t) = V Q {x,t) = V p (x,t-5t) 

$ Q (d,x;t) = -iV Q {x,t) = -iV p (x,t- St) 



(26) 



and we have 



${c,x x ;c,x 2 \t u t 2 ) = ^d ) x 1 ;d,x 2 ;ti,t 2 ) = ^[V Q (x 1 J 1 )V p (x 2) t 2 )+V p (x 1 j 1 )V Q (x 2 ,h)] (27) 
^{c,xi;d,x 2 ;t!,t 2 ) = ${d,x 2 ;c,x 1 ]t 2 ,t 1 ) = — [V Q (x 1 ,t 1 )V p (x 2 ,t 2 )~V p (x u t 1 )V Q {x 2 ,t^]- (28) 

One can see how the input amplitudes for the individual photons interfere when the two-photon wave function is written 
in terms of the input wave functions. As per the usual discussion, each photon is either reflected or transmitted by the 
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beam splitter giving four possible "outcomes" . For each outcome there is a corresponding wave function component in 
equations l|27l) - (|28[) . The $(c,x±;c,X2 ',ti,t2) and &(d,xi;d,X2',ti,t 2 ) components are the amplitudes for a detection 
event when both detectors are in the same output port. The amplitudes in (|28|) are those relevant to this experiment 
since we are interested in the coincidence detection rate where one photon is detected in each output port. The 
negative sign in (|28|l is the important feature that allows for destructive interference. A —i phase shift occurs when 
a photon wavefunction is reflected and the amplitude corresponding to a detection where both photons are reflected 
by the beam splitter accumulates a it phase shift relative to the amplitude of both photons being transmitted. The 
coincidence detection rate of the detectors is then proportional to 



\§{c,x\;d, x 2 ;ti,t 2 )\' 



\V p (x 2 ,t 2 -5t)\ + \V p {x 1 ,t 1 - St)\ \V p {x 2 ,t 2 )\ 
V p (a* , ii ) [V p {xx,ti- St)] * V p (x 2 , t 2 - St) [V p (x 2 ,h)Y 
[V p (x 1 ,h)]*V p (x 1 ,t 1 - St) [V p (x 2 ,t 2 - St)]*V p {x 2 M) . 



(29) 



The last two terms of Eq. 1|29[) are the exchange terms corresponding to the second two terms of Ijl9|) . As the P 
and Q (delayed) photon amplitudes overlap in the detection region, the exchange terms bring the detection rate to 
zero. At zero delay there is maximum interference and the coincidence detection rate is zero: The exchange terms 
completely cancel out the classical independent-particle amplitude for a two-photon detection. Similarity, evaluating 



\<b{c,X\; C,X2\tx,t2)\' t 



V p {xx,ti)\ \V p (x 2 ,t 2 - St)[ + \V p (xi,h - St)[ \V p (x2,h)\ 
+ V p (x 1 ,ti)[V p (x 1 ,ti~5t)Yv p (x 2 , h - St) [V p (X2,h)]* 
+ [V p (x 1 ,t 1 )]*V p (x 1 ,t 1 - St) [V p (x 2 ,t 2 ~ 8t)]*V p (x 2 ,t 2 ) 



(30) 



shows the constructive interference exchange effect that arises when measuring the coincidence detection rate of two 
detectors in the same output port. Assuming the photons are of finite temporal width, for large delay St we find 
|$(c, x\\ c, X2', ti, *2 ) | 2 = |$(c, x\, d, x 2 ; ti, t 2 )\ 2 . This may at first sight be surprising. According to the beam splitter 
transformations, each photon has a probability of 0.5 of being detected in port c, and a probability of 0.5 of being 
detected in port d. Since for large St we expect the photons to be independent, this would lead to a probability of 0.5 
that the photons are detected in separate exit ports, and a probability of only 0.25 that they would both be detected 
in port c, for example. However, it is important to recall that we have based the theory on the Glauber detection 
model where it is assumed that the interaction between the radiation field and detector is weak. When two detectors 
are located in the same port this allows for each photon in that port to interact with both detectors, in some sense 
double-counting the photons. This is the standard result of calculations based on the Glauber detection model. 

We now expand the discussion to include the non-separable case, where we allow for the possibility of "frequency 
entanglement". In the frequency entangled case Eq. I|18|) does not apply, and one cannot separately address the 
amplitudes of individual photons as previously discussed. Nonetheless, the photon wave function formalism can 
be used to describe frequency entangled two-photon interference in a slightly different way. Consider a two-port 
single-polarization gaussian entangled two-photon state as an input for the Hong-Ou-Mandel interferometer: 



\Pa,Qb) c = N(a) / dkdk'g p (k)g Q (k')e 



-at (k)al(k')\vac) 



N(a) 



2v^ 



die 



4 e 



-2iiUk„ t 



*(*)) 



where a is the entanglement width, 

l*(*)) = 

AT (a) = 



dkdk'g p {k)e iUkt g Q {k')e iuJ "' t al (k) a\ [k') \vac) , 

,2 , ^ ,2 2 (-fc+"V- 2 -fc„) 2 

dkdk' \g p (k)\ \g Q (k') e ^ 



(31) 
(32) 

(33) 
(34) 



and 2wfc is the peak entanglement angular frequency. One can think of this state as a superposition of separable 
two-photon states l)33|l . each with the two photons temporally displaced by a time t. In the superposition (|32|l the 
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amplitude of each component contain a phase factor that varies rapidly with this displacement time. The degree of 
the frequency entanglement determines the weighted distribution of these single-photon wave functions in time. It 
is in this form that one can discuss frequency entangled two-photon interference in terms of separable single-photon 
amplitudes. 

For the Hong-Ou-Mandel experiment with frequency entangled photons, one can describe the interference in terms 
of the separable input wave function amplitudes with a superposition of two-photon wave functions distributed in 
time. The two-photon detection amplitude becomes 



<&(c, xx;d, x 2 ; h,t 2 ) 



-4=^)4= 



dt [V p {xi,h - t)V Q (x 2 , h ~ t) - V Q (x!,ti - t)V p {x 2 ,t 2 - *)] 

(35) 



Comparing with (|28|l one can see that the phase shift of the reflected single-photon amplitudes once again allow for 
destructive interference. The detection rate is proportional to 



\<5>{c,x l ;d,x 2 -MM)\ 2 = T^l^)! 2 — 

61 7T 



-2iu ko (t-T) dtdT 



V p (x 1 ,t 1 - t)]*V p (x 1 ,t 1 - T) [V Q (x 2 ,t 2 - t)]*V Q (x 2 ,t 2 - T) 

+ [V p (x 2 , t 2 - t)\ * V p (x 2 ,t 2 - T) [V Q (x! M - t)] * V Q (x u ti - T) 

- [V p [ Xl ,ti - t)] * V Q (x 1 ,t 1 - T) [V Q (x 2 ,t 2 - t)]*V p (x 2 ,t 2 - T) 

- V p [x x ,h-T) [V* (x 1 ,t 1 -t)]*V c ^ (x 2 ,t 2 -T) [V p (x 2 , t 2 - *)] * 



(36) 



Here we have the possibility of exchange interference between amplitudes at different temporal displacements. Using 
this form it is possible to investigate the effect of frequency-entanglement on few-photon interference by examining 
the exchange terms. We will show that for a Hong-Ou-Mandel interferometer with a large enough delay such that the 
exchange terms are zero (no interference) with separable input photons, no amount of frequency-entanglement can 
introduce interference. We will begin with the assumption that the input wave functions can be well approximated as 
having finite temporal width, that is, if \a\ > (3 then V p (x, t)V p (x, t + a) = 0, where (3 is a photon width parameter. 
Given \St\ > ft, the exchange terms in (|29|l are zero and there is no interference in the separable case. To show 
that in this case there is no interference between temporally displaced two-photon wavefunctions, from examining the 
exchange terms in (|36|l we see we must show that for all (t, T) 



[V p (x 1 ,h+t)]*V p (xx ,h+6t + T) [V p (x 2 , t 2 + St + t)] * V p (x 2 , t 2 + T) = 0. 



(37) 



The proof is by contradiction. First we suppose that Eq. (|37|l is does not hold. With the finite temporal width 
assumption, this implies both \St + (T — t)\ < (3 and \6t — (T — t)\ < (3. But this is impossible to satisfy since \St\ > (3. 
Hence (13711 must hold and regardless of the entanglement, the exchange terms are zero if the initial relative photon 
delay, St, is larger than temporal width of the single-photon wave functions. One can see this in Fig. 2, which shows a 
pictorial representation of the exchange terms in (|36|l . In order for the exchange terms to contribute to the coincidence 
detection rate, all four single-photon wave functions must overlap at some value of (t, T). Regardless of the temporal 
displacements, the single photon amplitudes in the exchange terms do not overlap for St larger than temporal width of 
the photon wave functions. We will discuss later how this is not the case for the postponed compensation experiment, 
where the frequency-entanglement is more significant. 



B. Quantum eraser experiments 

The version of the quantum eraser we discuss here illustrates how the photon wave function formalism deals 
with multiple polarizations. The basic idea of the two-photon quantum eraser is as follows: It is possible to intro- 
duce distinguishing information in the input of an interferometer that would destroy second-order interference if the 
information were not "erased" at some stage before detection ^lj- I n our context, erasing distinguishing information 
is equivalent to producing overlapping wave functions in the detection region. 

Fig. 3 is a schematic diagram of a basic quantum eraser interferometer. Two orthogonally polarized photons in 
ports (a,b) interact with a beam splitter similar to the Hong-Ou-Mandel scenario. Since the photons are orthogonally 
polarized, it is possible to distinguish the final states of the two paths leading to a coincidence detection in ports c and 
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d. It is the presence of the 45 degree polarizing beam splitters in front of the detectors that causes the distinguishing 
information to be destroyed so that second-order interference can occur. The input state is described by 



\PaH,Q 



bV 



dkdk'g p (k)g Q (k')al H (k) a\ v (k') \vac 



(38) 



where, for example, a} aH (k) denotes a creation operator associated with a horizontally polarized photon of wave vector 
magnitude k in the propagation direction relevant for port a. As in the previous subsection, we define the input region 
wave functions. The two-photon wave function in the region just before the beam splitter has two components: 

$ HV (a,x 1 ;b,x 2 ;t 1 ,t 2 ) = ^ V H(b,x 2 ;a,x 1 ;t 2 ,t 1 ) = -^=V p (x 1 , h)V Q (x 2 , t 2 ), 

cf. <|23|) . Once again the relevant action of the beam splitter maps the input ports (a, b) to the output ports (c,d) in 
the following way: 



a bV ( k ) 



J.ui k A 



a dH ( fc ) - ia ln ( fc ) 
a[ v (k) - ia\ v (k) 



The components of the two-photon wavefunction just beyond the beam splitter are: 



(39) 



3>Hv(d,xi;d, xr,h,te) 
$ H v(c,x 1 ;c,X2;t 1 ,t2) 
$Hv(c,xi;d,x 2 ;ti,t2) 
<&VH{c,xi;d,x 2 ;ti,t 2 ) 



$vif(d, x 2 ; d, xv,te,ti) 
$VH{c,X2;c,x 1 ;t 2 ,t 1 ) ■ 
$VH(d,X2;c,xi;t 2 ,ti) 
$Hv(d,X2;c,xi;t2,ti) 



-L=V p (x l ,t 1 )V Q (x2,t 2 ) 
V p (x u t 1 )V Q (x 2 ,t2) 



2V2 



~^V p (x 1 ,ti)V Q (x 2 ,h) 

-j=V C >(x 1 ,tl)V P (x2,h), 



(40) 



cf. H25|l. To demonstrate that it is indeed the presence of the 45 degree polarizing beam splitters that change the 
system as to allow interference, we use the wave functions in this region to calculate the coincidence detection rate of 
the two detectors in absence of the polarizing beam splitters. In calculating this detection rate the square magnitude 
of <&Hv(c,xi;d,x 2 \ti,t 2 ) and $vh(c, x%; d, x 2 ; ti, t 2 ) are added separately, since 



,(2) 



(c,aji;d,x 2 ;*i,«a) = I \V p (x u t 1 )V Q (x 2 ,t 2 )\ 2 + i \V Q (x 1 ,t 1 )V p (x 2 ,t 2 )\ 2 . 



(41) 



No interference occurs if one removes the polarizing beam splitters since there are no exchange terms in (|41H . Now 
consider the experiment as shown. The effect of the 45 degree polarizing beam splitters is to map the ports (c,d) to 
the ports (e,f) as follows: 



a\ v (k) 
a dV ( k ) 



4A' 



V2 



\a\ + (k) + a\_ (fc)J , a\ H (fc) 
a) + (k) + a{_{k) ,a) dH {k) 



V2 



V2 



a\ + (k) - a\_ (k) 
4+ (k) - a{_ (fc) 



(42) 



where A and A' are the time intervals between the origins of the ports, and it is convenient to use delayed time 
coordinates t\( 2 ) — ti( 2 ) — A, i^) = ^1(2) — A'. Here we have used the + and — symbols to denote the +45 degree 
and —45 degree polarization bases respectively. The two-photon wave function component relevant to coincidence 
detection is 

$++(e,xi;f,X2;ti,t2) - $++(/, X2\e, x x ;t a ,*i) = — y= [V Q {xi 1 t 1 )V p {x 2 ,t 2 ) -V p {xiM)V Q {x 2 ,t 2 )] (43) 

which is the same form as the relevant wave function component in the detection region of the Hong-Ou-Mandel 
interferometer. Given that the photons initially have the same spectral decomposition and no relative delay, g p (k) = 
g®(k), and there is complete destructive interference and the coincidence detection rate of the detectors shown is zero 
as ® ++ (e,xi; f,x 2 ;ti,t 2 ) = 0. 
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C. The postponed compensation experiment 



The postponed compensation experiment was originally performed by Pittman el a/.|12|. The authors claim to 
have demonstrated an interference effect with two photons which do not arrive simultaneously at the beam splitter in 
an unbalanced Hong-Ou-Mandel type interferometer. They conclude that this effect can not be described in terms of 
the overlap of the individual photon wave packets on a beam splitter, and hence "two-photon interference" can not 
be considered as the "interference of two photons" . As we show below, the approach we have introduced here allows 
for a more precise formulation of such statements, thus both clearly identifying the physical insight they express and 
evaluating their validity. 

Fig. 4 shows the schematic diagram for the postponed compensation experiment. The input state consists of 
orthogonally polarized photons, with the vertically polarized photon in port b delayed by time t\ with respect to the 
horizontally polarized photon in port a. A relative delay, 72, is introduced in the horizontally polarized mode in the 
right arm of the interferometer. When T2 = 2ti the second delay compensates for the initial vertical photon delay in 
such a way as to create maximum interference for frequency entangled photons. While frequency entangled photons 
are used as the input state in this experiment, in order to understand this interferometer in the photon wave function 
picture we first suppose that one begins with the separable input state 



\P aH , P w ) = J J dkdk'g p (k)g p (k')al H (k) a\ v (*') \vac) 
for which the two-photon input wave function is given by 



(44) 



1 



§Hv[a,xi;b,X2\ti,t2) = $vH{b,x 2 ;a,xi;t 2 ,ti) = -j= [V p (xi, h)V p (x 2 , t 2 )\ 
The relevant action of the optics is to map the input ports in the following way 



(45) 



4» ( fc ) 

°lv (k) 



{ a t + (fc) - a \_ (fc)} e*»™ - i {al + (k) - a\_ (fc)} 
{al + (k) + at_ (fc)} e*"»" - i {a+ + (k) + a\_ (k)} 



which leads to the relevant two-photon wave function component in the region of the detectors 



3>++(c,xv,d,X2;ti,t2) = 



4\/2 



[V p (x 2 , t 2 - r 2 )V p {x 1 ,t 1 - ri) - V p (x 1 ,t 1 )V p (x 2 ,t 2 - ri)] , 



(46) 



(47) 



where t\{2) — £1(2) — A. When the interferometer is adjusted so that, t 2 = 2t\, one can see that the separable photon 
input yields no interference. The coincidence detection rate at maximum interference is proportional to 



\$ ++ (c,Xi;d,x 2 ;t 1 ,t2 



1 

32 



^(l 2 j2-27l)| IV^ixuh-Tt)] +\V P {x ll t 1 )\ \V P (x 2 ,t 2 -T 1 )\ 

- [V p (x 2 ,t 2 - 2T 1 )V p (x 1 ,t 1 - T 1 )]*V p {x 1 ,t 1 )V p (x 2 ,t2 - n) (48) 
-V p (x 2 ,t 2 - 2T 1 )V p (x 1 ,h - n) [V p (x 1 ,t 1 )V p (x 2l t 2 - n)] ! 

and the exchange terms are always zero, assuming the time delay t\ is much greater than the width of the wave 
functions. With the separable input state there are indeed no overlapping single photon amplitudes at the beam 
splitter, as Pittman et al. claim. However, one also does not observe any interference! 

It is only when the input photons are frequency entangled that one can measure an interference effect. Consider 
the frequency entangled input state 



P.,„.l\v) r = \-[<T) I I dkdk'g p (k)g p (k')e- { '" k+ " k y 



\ 2 

kg ) , , 

a an( k ) a lv (k')\vac) 



(49) 
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which leads to the relevant two-photon wave function component 



$ ++ (c,x 1 ;d,x 2 ;t 1 ,t 2 ) 



Km*)* 



-2iu)k t 



dt 



4\/2 ■ ' V 71 " 

The coincidence detection rate at maximum interference is then proportional to 



V P (X 2 , h -T2- t)V P {xi 2 h-Tl- t) 

-V p (x 1 ,h-t)V p (x 2 ,t 2 -n-t) 



(50) 



\$ ++ (c,x 1 ;d,x 2 ;t 1 ,t 2 ) 



> {t - T) dtdT 



[V p (x u ti - t)]*V p (x 1 ,t 1 - T) [V p {x 2 , t 2 - n - t)] * V p (x 2 ,t 2 -n-T) 

+ [V p {x 2 ,t 2 -2T 1 -t)]* V p {x 2 , t 2 - 2ri - T) [V p (x 1 , 1 1 - n - t)] * V p (x 1 , t r - n -{M) 

- [V p (x 1 ,h - t)]*V p {x 1 ,h -n-T) [V p {x 2 .t 2 - n - t)]*V p {x 2 ,t 2 - 2n - T) 

- V p (x 1 ,t 1 - T) [V p {xiM - n - t)]*v p (x 2 ,t 2 -n-T) [V p (x 2 ,t 2 - 2 n - 1)]' 

This expression is similar to that of the Hong-Ou-Mandel interferometer with frequency entangled photons (|36|l , but in 
this case the frequency-entanglement is necessary for interference. Given that the entanglement is sufficiently strong 
(1/(7 >> Ti), the exchange terms are non-zero, and one observes interference between the temporally distributed 
amplitudes of the two-photon wave function. A strong enough frequency entanglement introduces interference when 
there is no interference in the equivalent separable-photon input experiment. In order to show this we will again 
assume that the wave functions have a finite width: if \a\ > (3 then V p (x, t)V p (x, t + a) = 0, where /3 is the photon 
width parameter. Given |ti| > /3, the exchange terms in (|48|l are zero and there is no interference in the separable 
case. Examining the exchange terms in l|51|l we can observe interference in the frequency entangled case if 



[V P (x u h + t)] * V P (x 1 ,t 1 +n+T) [V p (x 2 , t 2 + Tl + t)] * V p {x 2 ,t 2 + 2ti +T) =/= 0. 



(52) 



With the finite width assumption Eq. (|52|l implies |ti + (T — t)\ < (3. Fig. 5 is a pictorial representation of how the 
temporally displaced single-photon wave functions overlap in the region of the detectors to produce non-zero exchange 
terms given \t\ + (T - t)\ < (3. 

In their discussion of this system, Pittman et al. state that the photons arrive at the beam splitter at much different 
times. Perhaps there exists some precise definition of what it means for "a photon to arrive at the beam splitter" for 
which this would be true. However, if one describes the frequency entangled photons according to the wave function 
picture discussed here, one finds that the wave function amplitudes of both the signal and idler are simultanously 
non-zero at the beam splitter. 

From a broader perspective, of course, we can agree with Pittman et al. that "the intuitively comforting notion of 
the photons overlapping at the beam splitter is not at the heart of the interference." One could - although Pittman 
et al. did not - perform an interference experiment where single-photon amplitudes do not overlap at a particular 
beam splitter (see Fig. 6). This would demonstrate, as Pittman et al. intended, that the notion of interference 
arising only when two single photons "meet" at a beam splitter is oversimplistic. Yet such an experiment would not 
demonstrate any limitation of a description of two-photon interference based on single-photon amplitudes. Even from 
this broader perspective we feel Pittman et al. go too far when they claim that "two-photon interference cannot be 
pictured as the interference between two single photons," at least insofar as it implies that there can be no general 
model of two-photon interference involving single-photon amplitudes. Indeed, we have presented such a model here. 
The coincidence detection rate does not explicitly depend on the amplitudes at a particular beam splitter or any 
other intermediate region, but rather the amplitudes in the detection region. We have shown in Eq. I|19|) that it is 
overlapping single-photon amplitudes in the detection region that gives rise to interference. 

Strekalov et q/.|l5| claim there are limitations of a single-photon wave packet approach for describing frequency 
entangled two-photon interference. Using the above theory to describe their experiment, one easily finds that the 
two-photon interference can in fact be understood in terms of temporally displaced pairs of single-photon amplitudes. 
In another paper, Kim et q/.|l6j| observe quantum interference between two temporally distinguishable pulses and 
discuss the limitation of a single-photon amplitude description. If one examines their experiment in terms of single- 
photon wave functions, one can see that the physics is the same as the postponed compensation experiment. The 
difference between the two experiments is that the temporal distribution is created by a series of pump pulses instead 
of a CW pump. The 50% visibility observed by Kim et al. is not surprising, at least in the photon wave function 
picture, as one can see it is simply postponed compensation interference with 50% probability. More recently, Kim 
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and Grice 17] report they have observed quantum interference in an experiment where the detected photons retain 
distinguishing information. Using the photon wavefunction theory to describe their experiment, one finds that the 
interfering single-photon amplitudes are overlapping in the region of the detectors. 

III. CONCLUSION 

Using a particular definition of the photon wave function, we have provided the formalism necessary for understand- 
ing second-order two-photon intereference in terms of individual photon amplitudes. The theory clarifies the idea of 
overlapping photons and replaces less rigorous explanations involving distinguishing information or Feynman paths. 
We have shown that the theory can be applied to both the separable and frequency entangled cases, and that in the 
latter it allows us to eliminate some of the confusion surrounding the interpretation of two-photon interference. This 
formalism shows how photon interference can be understood as an exchange effect by drawing an analogy to massive 
particle wavefunctions. For the systems discussed here, the presence of exchange terms in the coincidence detection 
rate expression is a necessary and sufficient condition for second order interference. Furthermore, if one considers only 
detectors with reasonable detection times and photon pairs with the same center frequency, second order interference 
is equivalent to the overlap of single photon wave functions in the detection region. 
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FIG. 1: A schematic diagram of the Hong-Ou-Mandel interferometer. Two like-polarized photons are incident on a beam 
splitter. The arrival of the photon in port b is delayed by 5t. Interference is observed in the coincidence detection rate of the 
two detectors. The origins of ports c and d are at the beam splitter, while the origins of ports a and 6 are at distance cA behind 
the beam splitter. 

T = t = 






FIG. 2: A pictorial representation of the Hong-Ou-Mandel exchange term photon wave functions in the region of the detectors. 
The dotted line indicates a two-photon wave function pair. In order for interference to occur, all four single photon amplitudes 
must overlap. Photon frequency entanglement allows for interference between temporally displaced two-photon wave functions. 
However, given that 5t is larger than the temporal width of the single-photon wave functions there is no interference between 
temporally displaced wave functions. 




FIG. 3: A schematic diagram of a quantum eraser. Orthogonally polarized photons are incident on a beam splitter. Interference 
in the coincidence detection rate is observed only in the region beyond the polarizing beam splitters PBS1 and PBS2. 




FIG. 4: A schematic diagram of the postponed compensation experiment. When ti — 2ti interference occurs in the coincidence 
detection rate despite a delay in mode 6. Frequency entangled photons are used in this experiment as no interference is observed 
with a separable photon state as the input. 
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FIG. 5: A pictorial representation postponed compensation exchange term photon wave functions in the region of the detectors. 
The dotted line indicates a two-photon wave function pair and P.C. denotes the postponed compensation apparatus. Once 
again, in order for interference to occur, all four single-photon amplitudes must overlap. No interference occurs in the separable 
case, but given a strong enough photon frequency entanglement the exchange terms contribute to the coincidence detection 
rate. One can observe interference in the frequency entangled case since the spatially displaced exchange term wave functions 
overlap if \r\ + (T — t)\ < f3, where /3 is the photon width. 




FIG. 6: A schematic diagram of an experiment designed to observe second order interference while the photons do not meet 
at the central beamsplitter. A separable photon input state is used. With n = T2 = r 3 maximum interference occurs in the 
coincidence detection rate. 



